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Abstract 

Based on Luttinger's formulation the complex optical conductivity tensor 
is calculated within the framework of the spin-polarized relativistic screened 
Korringa-Kohn-Rostoker method for layered systems by means of a contour in- 
tegration technique. For polar geometry and normal incidence ab-initio Kerr 
spectra of multilayer systems are then obtained by including via a 2 x 2 matrix 
technique all multiple reflections between layers and optical interferences in the 
layers. Applications to Co|Pt5 and Pts |Co|Pt5 on the top of a semi-infinite 
fcc-Pt(lll) bulk substrate show good qualitative agreement with the experi- 
mental spectra, but differ from those obtained by applying the commonly used 
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two-media approach. 
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I. INTRODUCTION 



Magneto-optical effects not only provide a powerful tool in probing the magnetic properties 
of solids, but are also of direct technological interest as phenomena to be used for high- 
density magneto-optical recording. Up to now, however, realistic theoretical investigations 
were lacking, because band-structure methods using supercells, [|] cannot provide an adequate 
description of layered systems, for which special computational techniques such as the spin- 
polarized relativistic screened Korringa-Kohn-Rostoker (SKKR) method have been designed. 
|^|-^ Furthermore, the absorptive parts of the optical conductivity tensor as obtained from the 
inter-band contributions, are not sufficient for magneto-optical Kerr spectra calculations, 
since also the dissipative parts have to be known. Hence, in supercell type calculations, besides 
the necessity to use the Kramers-Kronig relations, one has to include also the inter-band 
contributions by means of a semi-empirical Drude term. [1^ Only recently a better scheme was 
developed by two of the present authors , in which a contour integration is used to obtain the 



complex optical conductivity tensor as based Luttinger's formula, which in turn includes all 



inter- and intra-band contributions. [0 Combining this contour integration technique with the 



SKKR method, realistic inter- and intra-layer complex optical conductivities can be obtained 
for layered systems. 

Having evaluated the inter- and intra-layer optical conductivities, the magneto-optical Kerr 
spectra can then be calculated by using a macroscopical model such as, e.g., the two-media 
approach. []T1| Because the layered system contains more boundaries than just the interface 



between the vacuum and the surface layer, the two-media approach not fully includes the 
dynamics of the electromagnetic waves propagation in such systems. Since the pioneering work 



of Abeles in 1950, ||T^ several methods are known in the literature |]l6| , p!7| to treat multiple 



reflections and interferences using either a 2 x 2 matrix ||l8| , p!9| or 4 x 4 matrix ||20|- p^] technique. 
In the present paper the magneto-optical Kerr spectra of layered systems are evaluated for the 
most frequently used experimental set-up, namely polar geometry and normal incidence, by 
making use of the complex optical conductivity tensor and the 2x2 matrix technique. 

In Section |I| the theoretical background is reviewed briefly. Computational aspects are then 
summarized in Section In Section |IV| the two-media approach (Sect. |1V A| ) and the applied 



2x2 matrix technique (Sect. [IV C|) are viewed as two different macroscopic models of how to 
calculate magneto-optical Kerr spectra of layered systems. Particular emphasis has been put 
in Sect. [IV B| on the construction of layer-resolved permittivities in terms of the (macroscopic) 
material equation within linear response. This construction method combined with the 2x2 
matrix technique, allows one to determine self-consistently layer-resolved permittivities, see 
Sect. [IV (J4[ As an illustration ab-initio Kerr spectra of Co|Pt multilayer systems are presented 
and discussed in Section |V[ Finally, in Section |V| the main results are summarized. 



II. THEORETICAL FRAMEWORK 



A. Luttinger's formalism 



The frequency dependent complex optical conductivity tensor crioj) can be evaluated start- 
ing from the well-known Kubo formula and using a scalar potential description of the electric 



field. p3[ However, by using the equivalent ||2J] vector potential description of the electric field, 
one ends up with Luttinger's formula: |12[ 

S^^(u;) - S^j.(0) 



fku + iS 

with the current-current correlation function as given by 

/(^n) - /(em) 



V ^ /io; + zd + ( e„ - 6^) 



m,n. 



Here /(e) is the Fermi-Dirac distribution function, e^, e„ a pair of eigenvalues of the one- 
electron Hamiltonian, the are matrix elements of the electronic current operator (/i = 
X, y, z) and V the reference (crystalline) volume. 

The positive infinitesimal 5 implies the electromagnetic field to be turned on at t = —oo 



and hence describes the interaction of the system with its surrounding. |^ However, as can be 
seen from Eq. (^, b can be also viewed as a finite life-time broadening, which accounts for all 
scattering processes at a finite temperature. 

Luttinger's formula (|l|) and Eq. (§) have several advantages over the commonly used opti- 
cal conductivity tensor formula of Callaway. First of all, in contrast to the latter, Eq. ([l|) 
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simultaneously provides both, the absorptive and the dissipative parts of the optical conduc- 
tivity tensor. Hence there is no need for using the Kramers-Kronig relations in the Luttinger's 



formalism. On the other hand, as recently was shown |T3| Luttinger's formalism accounts for 
both, the inter- and the intra-band contribution on the same footing. Thus by using Eq. (|1|) 
in combination with Eq. (^, one does not need to include a phenomenological Drude term in 
order to simulate the intra-band contribution. |Q Furthermore, as was also demonstrated ||T^ 
Eqs. (|I]) and (Q) can be used for calculations in the static {u = 0) limit, provided the life-time 
broadening is kept finite (5 7^ 0). 



B. Contour integration technique 

Instead of evaluating the sums in the expression for the current-current correlation function 
(j^) over eigenvalues of the one-electron Hamiltonian, Ti^i,{uj) can be calculated by means of a 
contour integration in the complex energy plane. 

For the selection of a particular contour F, this technique |Tl] exploits the facts that with the 
exception of the Matsubara poles Zk = Sp + i {2k — 1) 5t {k = 0, ±1, ±2, . . . , and 6t = vt/cbT), 
P7| in both semi-planes the Fermi-Dirac distribution function of complex argument f{z) is 
analytical |25| and is a real function for complex energies z = e ± i6j situated in-between 

The latter property of f{z), e.g., is exploited by using 



two successive Matsubara poles. 



6j = 2Nj6T, where A^^i is the number of Matsubara poles included in F in the upper semi-plane 



and N2 in the lower semi-plane. |TT 



By applying the residue theorem, it has been shown [|TT| that equivalently to Eq. (Q) one 



has 



dz f{z) + C,^) 



dz f{z) j:^^{z - c, 



k=-N2 + l 



(3) 



such that 



d dz f{z)Y:^,y{z,z) -2i5T ^ Y:^,y{zk,Zk) 
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k=-N2 + l 



(4) 



where ( = hw + i6 and the kernel 

Z2) = Tr [J^ Giz,) r Giz2)] , (5) 

is related to the electronic Green function G{z). The auxiliary quantity Ti^^{zi, Z2) has already 
been used in residual resistivity calculations (u;, T = 0) of substitutionally disordered bulk sys- 
tems and magneto-transport calculations of inhomogeneously disordered layered systems. 
Pm Only recently, however, it has been shown, that Eqs. (3)-(|^) preserve all the advantages 
and features of Luttinger's formalism as was mentioned already above. 

In the present paper, S^^(zi,Z2) is evaluated in terms of relativistic current operators 
and the Green functions provided by the spin-polarized relativistic SKKR method for layered 



systems The optical conductivity tensor of a multilayer system is then given |31[ by 

N N 
p=l q=l 

with a^'^{(jj) refering to either the inter- (p 7^ q) or the intra-layer [p = q) contribution to the 
optical conductivity tensor. 

III. COMPUTATIONAL DETAILS 

In addition to the number of Matsubara poles considered, the optical conductivity tensor 
depends also on the number of complex energy points ra^ used in order to evaluate the energy 
integrals in Eqs. (3) and (^, on the number of k -points used to calculate the scattering path 
operators that define the Green functions and ± ^ + i5,z) for a given energy z. 

Recently, the present authors have proposed two schemes to control the accuracy of these z- 
and /c-integrations. 

The first of these schemes is meant to control the accuracy of the 2;-integrations along 
each contour part by comparing the results obtained from the Konrod quadrature 0,0, 
/C2n2+iS^jy(co'), with those from the Gauss integration rule, QnT,^y{uj). On a particular part 
of the contour T,^^{uj) is said to be converged, if the following convergence criterion: 

max lC2n,+l^i,y{'^) -Qn.^ixA'^) < ^ z ■, (7) 



is fulfilled for a given accuracy parameter Sz- 



The other scheme refers to the cumulative special points method, which permits to 
perform two-dimensional k -space integrations with arbitrary high precision. This method 
exploits the arbitrariness of the special points mesh origin. I^G] For a given (arbitrary high) 
accuracy the following convergence criterion has to apply 



max 



SnJ^^I.v{,z\ z) — Sn^_^^^y{z\ z) <e^, (8) 



for any complex energy z on the contour or Zk Matsubara pole. Here rii = 2*+^no (^o ^ N) is 
the number of divisions along each primitive translation vector in the two-dimensional A;-space 
and z' = z + (, z — (*. 

In the present paper, the optical conductivity tensor calculations were carried out for T = 
300 K, using a life-time broadening of 0.048 Ryd and N2 = 2 Matsubara poles in the lower 
semi-plane. Because the computation of cr^j,y(c<j) does not depend on the form of the contour. 



32 1 in the upper semi-plane we have accelerated the calculations by considering A^^i = 37 



Matsubara poles. The convergence criteria and (|^) were fulfilled for = ej: = 10 ^ a.u. 



IV. MAGNETO OPTICAL KERR EFFECT 

In the case of the polar magneto-optical Kerr effect (PMOKE), |]l4| the Kerr rotation angle 

= - ^ (A+ - A_) (9) 

and the Kerr ellipticity 

eK = - 10 

r+ + r_ 

are given in terms of the complex reflectivity of the right- (+) and left-handed (-) circularly 
polarized light 

f± = -^ = r± e^^± . (11) 

Here the complex amplitude of the reflected right- and left-handed circularly polarized light 
is denoted by S^'' and that of the incident light by £^^^; is the phase of the complex 



reflectivity r± and r± = \r±\. Eqs. (||) and ([10|) are exact, which can be easily deduced from 
simple geometrical arguments. However, in order to apply these relations, one needs to make 
use of a macroscopic model for the occurring reflectivities. 



A. Macroscopic model I: the two— media approach 

This simplest and most commonly used macroscopic model treats the multilayer system as 
a homogeneous, anisotropic, semi-infinite medium, such that the incident light is reflected only 
at the boundary between the vacuum and the surface (top) layer. 

In case of normal incidence the two-media approach provides an appropriate formula for 



the complex Kerr angle [|T4| 



^k = Ok- isk = , - , (12) 

which can be deduced from Eqs. (^) and (|T^ by assuming a small difference in the complex 
reflectivity of the right- and left-handed circularly polarized light. Because Im 5"xy (uj) usually 
is almost a hundred times smaller than Re (u), see Ref. [Q, the average complex refractive 
index of the right- and left-handed circularly polarized light is dominated by (Txx (^^) and hence 
the following direct formula results from Eq. (|T2|): 

4k « , ' . (13) 

1 - (Xxx [UJ) 

UJ 



with d-{uj) as given by Eq. (y). It should be noted that the "direct" formula in Eq. ([T3| ) 
was introduced by Reim and Schoenes [jl4[ in order to extract the optical conductivity tensor 



elements cXxx (uj) and o"xy (ui) from experimental PMOKE data. 



B. Macroscopic model II: layer— resolved permittivities 



Within linear response theory the Fourier transformed macroscopic material equations, 
38| averaged over the reference volume V, directly yields the total electric displacement 

^ [ d^rD{f,uj) = ^ [ d^r [ i{uj;r,f^)E{r',uj) , (14) 

V Jv V Jy Jy 



provided that the dielectric function e{uj;f, -F) and the Fourier components of the electric field 
E (r*, uj) are known. 

Using non-overlapping cells in configuration space (Atomic Sphere Approximation (ASA), 
applied in the present approach), the reference volume can be written as 



TV 



N 



V 



p=i 



p=i 



where N\\ is the number of atoms per layer (the same 2D lattice has to apply for each layer p), 
N the total number of layers and Qpi the volume of the ith atomic sphere in layer p. 

Assuming that plane waves propagate in a layer like in a two-dimensional unbound homo- 
geneous medium and that Dpi {f,uj) = Dp {f,uj), the integral on the left hand side of Eq. (|T^ ) 
can be written within the ASA as 



N 

d\ D (r, u) = N\\ XI 5Z 

p=l i 



1+6$: 



-1)'(A; + 1) [2 



k=l 



{2k + 3)\ 



'f^pSpi 

An 



2k 



(15) 



where P„ is the amplitude of the electric displacement, r?„ is the refraction vector. 



(16) 



Qp the wave vector (go = Stt/Aq refers to the propagation constant in vacuum) and Spi is the 
radius of the ith atomic sphere in layer p. Accordingly, the double integral on the right hand 
side of Eq. ([1^) reduces to 



N 



p,q=l 



V 



-1)' 



dr' (r') 
2n 



2k 



(17) 



^-(2/c + l)! VAo 

— * 

where £q is the amplitude of the electric field in layer q and e^^'^J (u;; r, r') is the dielectric 
function e (cu; r, F) at r G VLpi and P G Vtqj. 

In the case of visible light the wave vector dependence of the permittivity is negligeable. 
Therefore, after having substituted Eqs. (|I^) and (|T7p into Eq. (|T^), only the first term in 



the power series expansions has to be kept, which immediately leads to 



N 

E 



N 



where the inter- {p ^ q), intra-layer {p = q) permittivities are given by 



e^" (uj) 



(47r) I - ^^^2 I - ^^/)2 ^p^,gJ _ 



It should be noted that a similar result connecting the static current in layer p to the electric 
field in layer q is already known from electric transport theory in inhomogeneous or layered 
systems By using the relation Vp = eP (uj) Sp, the layer-resolved permittivities e^(a;) are 
then solutions of the following system of equations: 



N 



q=l 



p=l,...,N . 



By mapping the inter- and intra-layer contributions cr^'' (uj) to the microscopically exact optical 
conductivity tensor d-{uj), Eq. (^), onto the corresponding contributions of the permittivity 
tensor 



1 

N 



1 - (a;) 

UJ 



(19) 



one then can establish a well-defined macroscopical model for the evaluation of Kerr spectra. 



C. The 2x2 matrix technique 



1. Multiple reflections and optical interferences 



In contrast to the two-media approach the inclusion of all optical refiections and interfer- 
ences within a multilayer system assumes that each layer acts as a homogeneous, anisotropic 
medium between two boundaries and is characterized by a layer-resolved dielectric tensor 
{p=l,...,N). [0,0] 

As a first step the Fresnel or characteristic equation: p] | 



(/i, z/ = X, y, z) 



(20) 



has to be solved in order to determine the normal modes of the electromagnetic waves in a 
particular layer p. pT| Then by solving the Helmholtz equation for each normal modes, ET 







(/i,i^ = x,y,z) , 



(21) 



in 



the corresponding Spy components of the electric field in layer p are deduced. After having 



obtained the £pyS , the curl Maxwell equation [IS, IS 



(22) 



provides the amplitudes of the magnetic fields Tip for each normal mode in layer p. Here the 
Gaussian system of units has been used, Up is the refraction vector, as given by Eq. (|TB|), 
Up, which in an anisotropic medium is direction and frequency dependent. 



Finally, continuity of the tangential components of the electric and the magnetic field at the 
boundary between adjacent layers leads to a set of equations, which has to be solved recursively 
in order to determine the magneto-optical coefficients of the layered system, such as, e.g., the 
surface refiectivity. If no symmetry reduced quantities are used all the previous steps, Eqs. 



(|20|)-(22), have to be performed numerically by using - for example - the 2x2 matrix technique 
of Mansuripur. ||T8| , p!9 



Most frequently MOKE experiments are performed in polar geometry using normal inci- 
dence. Therefore in the following the 2x2 matrix technique of Mansuripur is confined to this 
particular experimental geometry. This reduction has the advantage that with the exception 
of the last step in which the surface refiectivity has to be evaluated, all the other steps can be 
carried out analytically. 

In the case of cubic, hexagonal or tetragonal systems and the orientation of the magnetiza- 
tion Mp pointing along the surface normal (z-direction), the layer-resolved permittivity tensor 
is given by 

(23) 



~p 



\ 



~p 

XX 


~p 

xy 





~p 

xy 


~p 

XX 











~p 

zz 



Assuming that e^^ — (p = 1, . . . ,N), the error introduced by this simplification - as easily 
can be shown - is proportional to the difference e^^ ~ which in turn is usually small enough 
to be neglected. If in polar geometry the incidence is normal. 



Tipx ^py 



for p=l,... ,N 
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the characteristic equation (^) provides four normal modes of electromagnetic waves in a layer 
p: 



Two of these four solutions are always situated in the lower half and the other two in the upper 
half of the complex plane. The first two solutions, h^pj and n^J , correspond to a "downward" 
(negative z -direction) and the other two, npj and fipz , sn "upward" (positive z-direction) 
propagation of the electromagnetic waves. These two different kinds of cases are given 

by 

ll'pz — V txx Itxy 

24 



- zexy 



and 



Tipz — \/ Cxx T ^Cxy 
f^pz — \/<;xx ^^xy 



(25) 



If in a given multilayer system a particular layer p is paramagnetic, its permittivity tensor 
is again of form shown in Eq. (p3|), with e^y = and e^^ — ^xx- ^^is case only two beams 
are propagating, namely those characterized by n^] = Up] = — v^HI and n^z = fi^z = 
Furthermore, since the vacuum is a homogeneous, isotropic, semi-infinite medium, in addition 

to 6xy = 0, exx — ^ZZ ~ 1- 

For each solution n^l^J (fc = 1, . . . , 4) of the characteristic equation (^Of), the electric field 
must satisfy the Helmholtz equation (|21|). Because not all of the equations are independent, 
these can be solved only for two components of the electric field keeping the third one arbitrary. 
Therefore, following the strategy proposed by Mansuripur, for beam 1 (n^J) and beam 3 {nj^^^ 



the corresponding Spx are chosen arbitrary, whereas for beam 2 (n^J) and beam 4 ( 

^pz ) , the ^py 

are arbitrary. For polar geometry and normal incidence, the solutions of the Helmholtz 

equation (|2T|) are given in Table | and the corresponding components of the magnetic field as 
obtained from Eq. (|2^) are listed in Table 
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2. Layer-resolved reflectivity matrix 



Numbering the layers starting from the first one on top of the substrate towards the surface, 
the surface layer has the layer index p = N , see Fig. |1|. The 2x2 reflectivity matrix TZp at the 
lower boundary Zp of layer p is given by: [0,0 



/ p(3) \ 



(3) 
px 

(4) 



7^„ 



( ~ . \ 







fL 



( c-a) \ 



•'px 
'(2) 



(26) 



see also the explicit discussion in the Appendix. The tangential components of the electric and 
magnetic field at a point just above the boundary then given by 



px 



y ^py j 

In ^ 



px 



y ^py j 



^ (J + 7^p) 



(2: - 7^p) 



/ ciX) \ 



y ^py J 
( ciX) \ 



(27) 



px 
(2) 



y ^py y 



where according to Tables H and p. 



1 I 



B. 



12 



V 



pz 



pz 



^(1) „-^(2) 



77, 



pz 



in- 



pz 



(28) 



and X is the 2x2 unit matrix. 

Using the lower boundary reference plane for the four beams in layer p — 1, the 

tangential components of the electric and magnetic field at a point just below the boundary 



Zp are of the form 



px 



y ^py y 



px 



y ^py y 



Bp'ii {Cp'ii - C^iiTZp-i) 



\ 



( cW 

'-'p-lx 

cm 
\ S-iy y 

c(2) 

\s-iy y 



(29) 



where 



C 



k,k+l 
p-1 



• -(fc+i) 



k = 1,3 



(30) 



with 



~(k) _ ~(k) , 



In here dp = Zp^i — Zp is the thickness of layer p, n^pli^ is defined in Eqs. ) and ( pSj) and go 



is the propagation constant in vacuum, see Sect. [IV B 



Based on Eqs. (^) and (p9D, the continuity of the tangential components of the electric and 
magnetic field on the boundary Zp implies that 



/ ^(1) \ / \ 



px 



cm 

( .(1) \ 

px 
(2) 



y ^py j 



K?12 /'/^12 _ ^-34 ^ 
\^p-\ ^p-\'^V-^j 



C-(2) 

yS-iy / 
/ 0(1) \ 

<-'p-lx 

\S-iy / 



such that by eliminating the electric field vectors, one immediately gets 



Vp_^{\^iip) = Bl' {x-n. 



■pj ) 



where 



T^p is therefore given in terms of 7^p_i by the following simple recursion relation: 
Hp = {Bf + Vp_,)-^ [Ef - Vp.,) p = 1, . . . , iV . 



(31) 



(32) 



In order to determine the refiectivity matrix TZ^ oi the surface layer, one has to evaluate all 
refiectivity matrices TZp for all layers below the surface layer. This requires to start the iterative 
procedure at the first layer (p = 1) on top of the substrate. But in order to calculate TZi, one 
needs to know the 2x2 matrix Vq corresponding to the substrate, see Eq. (|32|). This in turn, 
according to Eq. (^) is only the case, if the refiectivity matrix TZq of the substrate is available. 
In order to achieve this, one has to formulate the tangential components of the electric and 
magnetic field at z^ by taking into account that the substrate is a semi-infinite bulk without 
any boundaries and hence TZq = |T8|JT9|] Thus Vq = Bq^, which according to Eq. ( pHf ) requires 
to specify the permittivity of the substrate. 
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3. Surface reflectivity matrix 



In the vacuum region, since = 1, Cxy = 0, one has to deal with the superposition of only- 
two beams, namely that of the incident and reflected electromagnetic waves. These beams are 
related through the surface reflectivity matrix -Rsurf such that for polar geometry and normal 
incidence. 



<-'vac,x 



R 



■surf 



(i) 

vac.x 



vac,y 



^xx ^xy 



^xy ^xx 



(i) 

vac.x 



cii) 
y Ovac.y J 



(33) 



see the Appendix. Thus the tangential components of the electric and magnetic field at a point 
namely just above the interface between the vacuum and the surface, are given by 



where 



r / 

V 

v 

B. 



\ 



vac,x 



vac,y 



( cii) 
^ vac.x 



■surf J 



vac,y 



12 



/ 

v 



(-^vac ~^ -^vac-^surf) 



-N+l 

1 
-1 



y t.vac.y y 

/ cii) \ 
<-'vac,x 

cii) 
y <:'vac,y J 



(34) 



and B 



34 



V 







(35) 



According to Eqs. ( PU] ) and (|5^, the continuity of the tangential components of the electric 
and magnetic fields at the vacuum and surface layer interface, zn+i = 0, can be written as 



(T + -Rsurf) 



(i) 

vac,x 
(i) 

vac,y 



\ 



R12 _i_ 034 p 
vac ' -'-'var-'l'i 



34 



surf ) 



vac,x 



(i) 

vac,y 



( cii) \ 



0(2) 

/ .(1) \ 



E 



(1) 



ci-i) 



By eliminating from this system of equations the electric field vectors, it follows that 



J'n (1 + R : 



surf J 



r12 I p34 p 
^vac + ^..--R' 



surf 



where 



1-1 



(36) 



Thus for the surface reflectivity matrix one obtains 

-Rsurf = {^N - B'^^c) ^ [bIIc - ^n) ■ (37) 

The surface reflectivity matrix -Rsurf is therefore of the form given in Eq. (^3]), see also Appendix. 
In spherical coordinates, one immediately obtains the complex reflectivity of the right- and 
left-handed circularly polarized light as 

which in turn determines the Kerr rotation angle and ellipticity eK, see Eqs. (|]) and ([T0|). 



4- Self- consistent layer-resolved permittivities 



In order to calculate for a homogeneous, anisotropic layer p the corresponding dielectric 
tensor (pSj) from the inter- and intra-layer permittivities defined in Eq. (|l^), the following 
linear system of equations 



v 



IP IP 

XX xy 



-xy -^xx 



px 



y ^py I 



N 

E 

5=1 



v 



~pq ~pq 
XX xy 



7pq 2Pq 
xy XX 



qx 



has to be solved, see Eq. (pISl). Here for one can take the following Ansatz: 



px 



y ^py j 



where due to Eq. ( pO| 



px 



A 



y ^py j 



( ciX) \ 



■'px 
'(2) 



y ^py j 



( _. -(k)dp 



V 







. -(fc+l)dp 



k = l,3 . 



By using the continuity equation of the tangential components of the electric field at the 
boundaries, see Eqs. (pT]), ( pQ]) and (0), one then obtains the layer-resolved permittivities as 
a weighted sum of the inter- and intra-layer permittivities defined in Eq. (|T9|) : 



"'XX xy 



_~p ~p 

xy XX 



(3J 



q=l 



where 



-1 



w^p. = ^ I n ^^+'^ 1 1 n ^v+k 



(39) 



,fe=0 



with 



Wp+fc = (J + 7^p+fc) (Cp^2 ^ + clXjjZp+k) \ k = i,...,N-p 



and 



(Cf)U(0^7^J(C- + Cf7^,)-^ , A; = 0. 



Because the 2x2 matrices Wp+fc contain T^p, C^^ and C^^, which in turn depend on layer-resolved 
permittivities e^j, (cj), Eq. (^) has to be solved iteratively. 

The self-consistent procedure can be started by putting all 2 x 2 weighting matrices Wpq in 
Eq. (|39D to unity, i.e., by neglecting the phase differences of the electromagnetic waves between 
the lower and upper boundaries in each layer p, 



N 

E 

q=l 



(40) 



These quantities e^j^ ('^)*'°^ can be used to calculate TZf"* in terms of Eqs. (|31|) and (|32|). Improved 
layer-resolved permittivities follow then from Eq. (pH]). This iterative procedure has to be 
repeated until the difference in the old and new layer-resolved permittivity of layer p is below 



a numerical threshold e. 



max 



< e 



p ■ 



(41) 



V. RESULTS AND DISCUSSIONS 



/^From experiments is known that Pt substrates prefer an fcc(lll) orientation. [^] There- 
fore in the present contribution, the calculations for the layered systems Co|Pt5 |Pt(lll) and 



Pts |Co|Pt5 |Pt(lll) have been performed, with five Pt layers serving as buffer |^3| to bulk fee 
Pt. 

The Fermi level in Eqs. (3) and (H) is that of paramagnetic fee Pt bulk (lattice parameter 



of 7.4137 a.u.), which also serves as parent lattice, H3] i.e., no layer relaxation is considered 
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A. Paramagnetic fcc(lll) Pt substrate 



As mentioned above in order to determine the surface reflectivity the permittivity tensor 
of the semi-infinite substrate has to be evaluated. As can be seen from Fig. 0, the xx-element 
of the permittivity tensor of the fcc(lll) Pt substrate shows a rather simple photon energy 
dependence. The real part of the permittivity exx(i^) has a peak around 1 eV, while the 
imaginary part of exx(t^) exhibits an almost perfect hyperbolic frequency dependence. The 
strong decay of exx(ci;) for photon energies in the vicinity of the static limit {uj = 0) can be 
easily understood in terms of the Eqs. ( ]TP| ) and (^0]), see also Ref. : for c<j the real part 
of ix^{uj) must tend to minus infinity whereas the imaginary part has to decrease. 

The xy-element of the permittivity tensor for fcc-Pt(lll) is identical zero over the whole 
range of optical frequencies, a functional behavior that of course does not need to be illustrated. 



B. Self consistent layer resolved permittivities 

In terms of the substrate and the zeroth order layer-resolved permittivities, see Eq. (^0[), 
the iterative determination of the surface refiectivity matrix described above, provides also 
self-consistent, layer-resolved permittivities i^xi^) ^ very efficient manner: in less than five 
iterations an accuracy of Sp = 10^^^ for each layer p, see Eq. (|^), can be achieved. 

In order to illustrate this procedure, in Fig. |^ the imaginary part of the relative difference 
between the self-consistent and zeroth order layer-resolved xx-element of the permittivity 
tensor for Co|Pt5 |Pt(lll) with and without Pt cap layers is displayed. 

This relative difference is to be viewed as the relative error made by using according to 
Eq. (^) the 2x2 matrix technique with zeroth order layer-resolved permittivities. As can 
be seen from Fig. 0, this relative error is layer, frequency and system dependent. The higher 
the photon energy and the bigger the layered system the less exact are the zeroth order layer- 
resolved permittivities. However, for relatively small layered systems, the relative error made 
by using only zeroth order permittivities is typically below 5 % for exx(i^) and less than 20 % 
for exy(co'). However, the resulting relative error in the Kerr rotation angle and the ellipticity as 
calculated by comparing the spectra corresponding to the self-consistent and to the zeroth order 



18 



layer-resolved permittivities, is always less than 1 %. Therefore, Eq. pO|) can be considered as 
reasonably good approximation for layer-resolved permittivities. 



C. Polar Kerr effect for normal incidence 

The systems investigated in here refer to a Co mono-layer on the top of a fcc-Pt(lll) 
substrate, considering also the case of three Pt cap layers. As already mentioned, five Ft layers 
serve as buffer to the semi-infinite host in order to ensure that the induced magnetic moments 
decrease monotonously to zero into the paramagnetic Pt substrate. 

The ab-initio Kerr spectra obtained from self-consistent layer-resolved permittivities by 
applying the 2x2 matrix technique are shown in Fig. ^. Usually, in experiments Pt cap 
layers are deposited on top of Co in order to prevent the oxidation of the surface. [H4| By 



performing a separate, magnetic anisotropy calculation, we have found that Co|Pt(lll) 
exhibits perpendicular magnetization only in the presence of Pt cap layers. Therefore, for the 
polar Kerr spectra of the CojPts |Pt(lll) system shown in Fig. ^ the polar geometry, namely 
the perpendicular orientation of the magnetization is imposed. 

Analyzing in Fig. ^ the Kerr spectra of the capped and that of the uncapped systems, several 
differences can be observed. The negative peak in the Kerr rotation angle at 3 eV in the 
spectrum of CojPts |Pt(lll) almost disappears from the spectrum in the case of the capped 
layered system. In the Kerr ellipticity the zero location at 2.8 eV observed for the uncapped 
system is shifted to 2.5 eV for the capped system system and simultaneously the infrared (IR) 
positive peak is shrunk and moved towards lower photon energies. Besides these features, the 
sign of the ultraviolet (UV) peak in both, the Kerr rotation and ellipticity spectra is changed, 
when the Co surface layer is capped. This particular feature can also be observed in the Kerr 
spectra obtained by using the two-media approach, see Fig. |^. In the two-media Kerr rotation 
spectrum, the negative IR peak at 2 eV in Co|Pt5 |Pt(lll) is shifted to about 1 eV in case of 
the capped system. 

Comparing the spectra in Fig. ^ with those in Fig. ^ it is evident that the theoretical 
Kerr spectra depend indeed on the macroscopic model used to describe the propagation of 
electromagnetic waves in the system. Because the systems investigated in here are much smaller 

IP 



than those used in experiments, p^ , ^ a strict quantitative comparison with experimental data 
cannot be made. However, a quahtative comparison based on the well-known, general features 
of the Co|Pt experimental Kerr spectra is still possible: ||4^ the Kerr rotation angle shows (a) a 
small negative IR peak at 1.5 eV, which decreases in amplitude with decreasing Co thickness, 
and (b) a high and broad negative UV peak, which moves from 4.1 eV to 3.9 eV for increasing 
Co thickness. The Kerr ellipticity is characterized by (a) a shift of the zero location at 1.5 eV 
(pure Co film) to 3.7 eV with decreasing Co thickness, (b) a positive peak around 3 eV and (c) 
a shift of the minimum at 4.9 eV in pure Co film towards higher photon energies. 

In the Kerr rotation spectra of the capped system shown in Fig. ^, there is no negative IR 
peak around 1.5 eV, but a negative UV exists at 5.5 eV. The Kerr ellipticity spectra of the 
capped system in Fig. ^ has a zero location at 2.5 eV and positive peaks show up at 0.5, 3.5, 
4 and 5 eV. These features suggest that in the case of Pta |Co|Pt5 |Pt(lll) the Kerr spectra 
obtained by applying the 2x2 matrix technique are typical for Co|Pt layered systems. 

A similar investigation of the Kerr rotation spectra in Fig. |^, reveals that for the capped 
system there are two negative IR peaks at 1, 1.5 eV and a negative UV peak around 5 eV. The 
Kerr ellipticity for the capped system in Fig. ^ shows a zero location at 1 eV (1.5 eV in case 
of a pure Co film), two positive peaks at 4 and 5 eV and a small negative peak around 3 eV. 
All these features make the Kerr spectra of Pts |Co|Pt5 |Pt(lll) described via the two-media 
approach to resemble those of a pure Co film rather than those of a Co|Pt layered system. 

Previous results obtained by applying the 2x2 matrix technique using as substrate per- 
mittivity that of the last Pt-layer below the Co one, have shown similar characteristics in 
the Kerr spectra. Hence, these features cannot be ascribed to the presence of the substrate. 



In an other contribution, it was shown that the optical conductivity of these systems is 
dominated by the contributions arising from the polarized Pt-layers. Therefore, the pure Co 
film-like spectra obtained for Pt-^ |Co|Pt5 |Pt(lll) within the two-media approach can be seen 
as an indication that a layered system cannot be approximated by a homogeneous medium in 
which no optical reflections or interferences occur. 
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VI. SUMMARY 



We have used the 2x2 matrix technique for the most frequently used experimental set-up, 
namely for polar geometry and normal incidence. This technique allows one to account for 
all multiple reflections and optical interferences in a semi-infinite layered system. The Kerr 
rotation angle and ellipticity can be directly obtained from the iteratively calculated surface 
reflectivity matrix, which in turn can be used to determine layer-resolved permittivities self- 
consistently. For free surface of layered systems realistic ab-initio Kerr spectra are obtained 
using the inter- and intra-layer conductivities as given by Luttinger's formula within the spin- 
polarized relativistic screened Korringa-Kohn-Rostoker method. 

A comparison of the theoretical Kerr spectra of Co|Pt5 |Pt(lll) and Pta |Co|Pt5 |Pt(lll) as 
obtained by applying the 2x2 matrix technique and the two-media approach indicates that 
the former technique provides typical results for layered systems, whereas the later approach 
tends to generate spectra specific for homogeneous films on top of a substrate. 



This work was supported by the Austrian Ministry of Science (Contract No. 45.451), by 
the Hungarian National Science Foundation (Contract No. OTKA T030240 and T029813) and 
partially by the RTN network "Computational Magnetoelectronics" (Contract No. HPRN- 
CT-2000-00143). 

APPENDIX: SYMMETRY OF REFLECTIVITY MATRICES 

Since for a semi-infinite substrate, TZq = 0, Vq = B^"^, with B^'^ as given by Eq. (^81) , 
according to Eq. (|3^) , the reflectivity matrix of the first layer on top of the substrate is given 
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ri 



~(i) 



n 



Iz 



Oz 



;^(2) ^^(2) 



n 



Iz 



Oz 



~(2) , ~(2) 
'^Iz + ^Oz 

Assuming that all p — 1 reflectivity matrices are of this diagonal form, namely 



T, 



j = i,...,(p-i) 



by taking into account Eqs. (^) and (|30D , Eq. (|3l|) immediately yields 



-id, 



where 



C?p_l 



.-{1) .-(3) 
"■p-lz ^.,(1) _i,-(3) 



;;(2) 



• -(4) 
"P 1 . 



n, 



p-i 



p— Iz 



(2) 



p-1 _|_ g ^fp-lf 



p-1 



The reflectivity matrix of layer p as obtained from the recursion relation ( p2D is found to be 
also diagonal 



rp 



where 



' t-pz Wp_ I 
fT'pz ~\~ dp—1 
pz fJ-p-l 



'tpz "p_l 

i.e., all the layer-resolved reflectivity matrices TZj { j = 1, . . . , N) are diagonal matrices as was 
anticipated in Eq. (]26|). 

In terms of the diagonal reflectivity matrix of the surface layer T^^r, and A and B]^ as given 
by Eq. (|8D, Eq. ( H) reduces to 



1 / ifN -In^ 



N 



J'n */jv 
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where 



Jn — (I'n ~ d-N 



/jy — c?^ + d 



N 



By using jF/v together with the matrices defined in Eq. (^) in Eq. (|37|) , the resuhing surface 
reflectivity matrix is of the form anticipated in Eq. (|33D, i.e. 



surf 



^xx ^xy 



' ^xy ^xx 



where 



N 



-Ai- 



N 
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TABLES 



k 


1 


2 


3 


4 


Opx 


arbitrary 


LCpy 


arbitrary 


ICpy 


Cpy 


lOpx 


arbitrary 


iOpy^ 


arbitrary 


cik) 
Cpz 















TABLE L Solutions of the Helmholtz equation (|2^) for polar geometry and normal incidence ne- 
glecting the difference in the diagonal elements of the layer-resolved permittivity. Ej^l is the amplitude 
of the electric field in layer p for beam k. 



k 


1 


2 


CO 


4 


-yjik) 
1 i-px 




py V 


iEpx ^ 6xx ~l~ ^^xy 


_c(4) /~p _ .-rp 
py V XX ^xy 


rxpy 


C'px Y txx \ 'txy 


't-py Y txx 'txy 


C-px Y txx \ 'fcxy 


•c-(4) /~p .-rP 
tt-py Y txx 'txy 


■yj(k) 
1 i-pz 















TABLE IL Solutions of the curl Maxwell equation ( p2|) for polar geometry and normal incidence 

{k\ 

neglecting the difference in the diagonal elements of the layer-resolved permittivity e^. Upp, is the 
amplitude of the magnetic field in layer p for beam k. 
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FIGURES 
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FIG. 1. The macroscopic model used for a layered system within the 2x2 matrix technique for 
polar geometry and normal incidence. The x axis is perpendicular to the plane of the figure, g^*) is 
the incident and g is the reflected wave vector. M denotes the total spontaneous magnetization of 
the system. 
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FIG. 2. The permittivity for fcc(lll)-Pt bulk as a function of the photon energy u. The real 



part of the permittivity is denoted by full circles, the imaginary part by open circles. 
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FIG. 3. Imaginary part of the relative difference between the self-consistent and zeroth order 
layer-resolved xx-element of the permittivity tensor as a function of the photon energy u for fee 
Co|Pt5|Pt(lll) (top) and Pta |Co|Pt5 |Pt(lll) (bottom). The data represented by full (open) circles 
correspond to the first, squares to the second and diamonds to the third Pt layer on top of (under) 
the Co layer (stars). Open triangles down (up) denote the first (second) Pt layer data on top of a 
paramagnetic fcc-Pt(lll) substrate. 
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FIG. 4. The magneto-optical Kerr rotation angle %(a;) and ellipticity eK(a;) for polar geometry 
and normal incidence as a function of the photon energy uj obtained by applying the 2x2 matrix 

technique for the self-consistent layer-resolved permittivities of fee Co|Pt5 |Pt(lll) (open circles) and 
Pt3|Co|Pt5|Pt(lll) (full circles). 
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FIG. 5. As in Fig. ^, but here the Kerr spectra was obtained by applying the two-media approach 
for fee Co|Pt5|Pt(lll) (open circles) and Ptg |Co|Pt5 |Pt(lll) (full circles). 



